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ABSTRACT 

We analyze 1298 merging galaxies with redshifts up to 2; = 0.7 from the Canada-France-Hawaii 
Telescope Legacy Survey, taken from the catalog presented in Bridge et al. (2010). By analyzing the 
internal colors of these systems, we show that so-called wet and dry mergers evolve in different senses, 
and quantify the space densities of these systems. The local space density of wet mergers is essentially 
identical to the local space density of dry mergers. The evolution in the total merger rate is modest 
out to z ~ 0.7, although the wet and dry populations have different evolutionary trends. At higher 
redshifts dry mergers make a smaller contribution to the total merging galaxy population, but this 
is offset by a roughly equivalent increase in the contribution from wet mergers. By comparing the 
mass density function of early-type galaxies to the corresponding mass density function for merging 
systems, we show that not all the major mergers with the highest masses (MstcUar > lO-'^^M©) will end 
up with the most massive early-type galaxies, unless the merging timescale is dramatically longer than 
that usually assumed. On the other hand, the usually-assumed merging timescale of ^ 0.5 — 1 Gyr 
is quite consistent with the data if we suppose that only less massive early-type galaxies form via 
mergers. Since low-intermediate mass ellipticals are 10-100 times more common than their most 
massive counterparts, the hierarchical explanation for the origin of early- type galaxies may be correct 
for the vast majority of early- types, even if incorrect for the most massive ones. 



1. INTRODUCTION 

Hierarchical structure formation models suggest that 
galaxy mergers play an important role in galaxy mass as- 
sembly, but quantifying that role has remained challeng- 
ing. The galaxy merger rate is generally parametrized 
by a power law of the form (1 -I- z)™, and the value 
of the exponent of this parametric form has been used 
to place constraints on how much mass is assembled 
via major galaxy mergers. Large variations in m 
are fou nd in the l iterature, ranging from m ~ to 
4 (IBundv et al. 2004, 2009; Consclicc e t al 



Guo fc \Vh ite 2008 ; iLin et a l. 2004; Bridge ctraH 
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2007 



Lin et all 12008: L otz et all 120 08b: Jogcc et al. 200l). A 



recent study by [Bridge et al.l ((20101) analyzed these pub- 
lished merger rates and concluded that, overall, there is 
a general agreement that the merger rate at intermedi- 
ate redshifts (0.2 < z < 1.2) does evolve, although the 
constraints on m remain fairly mild. Bridge et al. (2010) 
rule out m < 1.5 (i.e. flat or mild evolution) and sug- 
gest that the wide range of m reported in the literature 
is due to a combination of factors, including variation in 
the redshift ranges being probed, small sample sizes in 
some of the surveys, and cosmic variance. 
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The traditional approach used to explore the merger 
history of galaxies has been to estimate the fraction of 
merging galaxies relative to the total galaxy population. 
This approach has the benefit of simplicity, but it is 
arguable that a more physically interesting quantity is 
the evolving space density of mergers, rather than the 
merger fraction. Space densities are absolute measure- 
ments rather than relative measurements, and in that 
sense can stand on their own. Furthermore, space den- 
sities can be corrected for luminosity biases and other 
sources of incompleteness in a straightfor ward manner 
by using the standard 1/T4iax formalism (|Feltenl [19771 : 
ISchmidtl 119681) ■ Therefore our main aim in the present 
paper is to chart the evolving space density of merg- 
ers. Simil ar work has b een u ndertaken by Ljn et aLl 
(|2008D and IBundv et al.l ()2009[) . ahhough these papers 
used pair counts to select merging galaxies, while our 
approach is based on morphologica l selection . Our anal- 
ysis is thus quite co mplementary to ILin et"aLl (|2008[ ) and 
IBundv et all (|2009f). 

An important subsidiary goal of the present paper is 
to chart the differential merging history of color-selected 
sub-classes of merging galaxies. In recent years a host 
of observations have shown the evolutionary histories of 
galaxies in the so-c alled 'red sequence' and and 'blue 
cloud' are different ((Bundv et al.l 120091: Ide Ravel et al.l 



[200l ILin et alJ [20081: IWillmer et all [20061) . This has 
led to the idea that it is important to distinguish be- 
tween mergers that result in significant star-formation 
('wet mergers') and those which merely re-organize ex- 
isting stellar populations ('dry mergers'). Wet mergers 
are typically ass ociated with ga s-rich systems that trigger 
star formation (iBarton et al.l [20 00: Bridge et al. 2007; 
Lin et all [20071: lOverzier et al.ri200& .Pi Matteo et al.l 



20081: IRodrig uez Zaurm et al. r i2010D . cause quasar ac- 
tivity (Hopkins et al.' '2006") and turn disk galaxies into 
elliptical galaxies (Toomro & Toomrc 1972). On the 
other hand, dry mergers ( Bell et al.l 120061 : Ivan DokkumI 
[2005h may be responsible for the assembly of massive 
(Msteiiar ^ 10^^ M©) red galaxies which are observed 
in surprisingly high abundance at z 1. However, the 
importance of these wet and dry mergers in the forma- 
tion of red sequence galaxies is still not clear (iBell et al 
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Brown et al.ll20q7l: iBundv et al.ll2009l: iFaber et al. 
Lin et al.l 120081: IScarlata et al.ll200'7[) . and it is of 



interest to determine if wet and dry merging systems ex- 
hibit similar evolutionary trends as a function of cosmic 
epoch. 

A plan for this paper follows. In Section 2 we describe 
the observations, galaxy properties and merger identi- 
fications for our sample. This section is essentially a 
brief recapitulation, presented for the convenience of the 
reader, of the comprehensive description of the data set 
given in Bridge et al. (2010). In Section 3 we define 
our methodology for defining 'wet' and 'dry' mergers on 
the basis of both resolved colors and integrated colors. 
Section 4 is the heart of the paper, where we present the 
merger fraction, number density and stellar mass density 
for dry and wet mergers. Since recent work shows that 
the formation rate of massive elliptical galaxi es through 
dry mergers is dependent on stellar masses (|Bell et al.l 
[20071: IBundv et al.l[20Ql IKhochfar fc"Siig [200^ 7 in this 
section we also investigate the evolving space density of 
mass-segregated samples of wet and dry mergers. Our 
results are discussed in Section 5, and our conclusions 
presented in Section 6 

Throughout this paper, we adopt a concordance cos- 
mology with Ho=70 km s~^ Mpc~^, f^M = 0.3, and 
fl\ = 0.7. All photometric magnitudes are given in the 
AB system. 

2. OBSERVATIONS 

As has already been noted, a detailed description of 
the selection strategy for (and basic properties of) the 
galaxies analyzed in the present paper has already been 
presented in Bridge et al. (2010). The reader is referred 
to that paper for details beyond the outline presented 
here. 

2.1. Data 

The data in this paper come from two of the Canada- 
France-Hawaii Telescope Legacy Survey (CFHTLS) deep 
survey fields. These fields (denoted Dl and D2) to- 
gether cover an area of 2 square degrees. The CFHTLS 
deep survey has high-quality broad-band photometry in 
five bands (u*, g', r',i',z') and the depth of the survey 
ranges from 26.0 {z') to 27.8 {g'). The optical images 
used to derive galaxy properties and morphological clas- 
sifications were stacked by the Elixir image processing 
pipeline (|Magnier k. Cuillandre ,2004 ) to produce deep 
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Fig. 1. — A montage of i'-band image of six typical galaxy merg- 
ers. Each stamp has the dimension of 100 X 100 kpc. The catalog 
name and the photometric redshift of the merger is shown in the 
bottom left corner of each stamp. 

optical stacks with precise astrometric solutions. The 
typical seeing for the final stacks is 0.7"-0.8" in the i'- 
band. The source extraction and pho tometry were per- 
forme d on each field using SExtractor (jBertin fc Arnouta 
UME) in dual image mode. The source detection was per- 
formed in the i' filter {i' ^ 26.3). A bad pixel mask was 
applied to the image prior to running the program to 
avoid noisy or contaminated regions caused by spikes or 
halos of bright stars. The total area masked is less than 
10% for each field. 

2.2. Merger identifications 

Merging galaxies were selected visually, with multiple 
cross-checks on the visual classifications, and using sim- 
ulations to characterize detection thresholds for features 
that are signatures of mergers. Visual classifications are 
labor intensive and somewhat subjective, but remain the 
best method presently available if accuracy is the ulti- 
mate goal. Interacting galaxies are defined as systems 
with a tidal tail or bridge. All galaxies down to an 
*«ega — ^2.9 mag (~ 27,000) were inspected resulting 
in a final sample of 1298 merging galaxies. The merger 
identification rate for galaxies with i'^^g^^ < 21.9 mag 
is estimated to be > 90%, and drops rather quickly 
at fainter magnitudes. Therefore in the present paper 
the iygga = 21.9 mag is used as the limit for merger 
identification. The redshift completeness limit was es- 
timated by artificially redshifting low-redshift galaxies. 
For this completeness test, a group of galaxy mergers 
with redshift ranges from z = 0.3 to z = 0.45 and 
Mg < —21.0 mag were selected. They were artificially 
redshifted to higher redshifts after accounting for the k- 
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Fig. 2.— Colors of ~ 2200 field galaxies in the CFHTLS survey. 
The red filled circles and blue triangles represent the visually classi- 
fied elliptical galaxies and spiral galaxies, respectively. The dashed 
line with a slope of 0.13 is derived from the red sequence fitting 
of 115 visually-classified elliptical galaxies with spectroscopic red- 
shifts from the COSMOS; two equally spaced parallel dotted lines 
are used for merger classification test (described in the text). The 
horizontal dashed line represents rest frame color g' —i' = 1.29 mag, 
and is also used for the same test. The green circles and triangles 
represent elliptical galaxies and spiral galaxies with MIPS 24/im 
detections. As described in the text, we assign the MlPS-detected 
galaxies with g' — i' color greater than the color bimodality line 
to the blue population. The results of merger classification test is 
shown in Figure |3] see text for more information. 

correction, change in angular size and surface brightness 
dimming. After this step the merger identification was 
conducted once again and found to be acceptably high 
out to z ~ 0.7, where ^ 85% of the redshifted merg- 
ers could be still be identified as mergers, and the false 
positive contamination fraction remains minimal. This 
redshift limit is particularly important because we uti- 
lize a l/Vmax weighting to correct the Malmquist bias 
and to compute number densities, li is important to note 
that Vniax of mergers presented in this paper is defined as 
the maximum volume over which mergers can be iden- 
tified as such, and not the maximum redshifi, at which 
a given galaxy's integrated magnitude remains above the 
detection threshold. As described in the next section, we 
calculated the Knax value from the Zmax provided by the 
Z-Peg code which denotes the maximum redshift that 
the template SEDs is fainter than the observed limiting 
magnitude (i^ega ^ 21.9 mag). 

2.3. Galaxy properties 

The CFHTLS survey has high quality five broad band 
photometry which makes the derivation of accurate pho- 
tometric redshifts, ages and stellar masses possible. The 
galaxy properties were derived by comparing the spec- 
tral energy distributions (SEDs) obtained from observed 
fluxes to a set of template SEDs. The best-fit SEDs 
were determined through a standard minimum fit" 
ting between the template SEDs and the observed fluxes. 
The template SEDs were co mputed by the PEG ASE-II 
galaxy evolution co de (Fioc & Rocca-Volmcrangc 1997 



Le Borgne fc Rocca- Volmcrangc 2002; Lc Borgnc ct alJ 



20041 ) and were integrated through the CFHT filters. 



The SED fits were undertaken using the Z- Peg code 
(iBolzonella et all 120001 : iLe Borgne &: Rocca-Vo lmcrangj 
I2002D and details are described in [Bridge et al.l ((201^ 
The photometric accuracy is determined by compar- 



ing the derived photometric redshi fts to the spectro- 
scopic redshifts in th e SNLS sample (|Howell et al.l 120051 : 
iBronder et all [20081 ) . The accuracy of the photometric 
redshift down to i ~ 22.5 is ctaJU + z) = 0.04. The 
stellar mass for each merger was also estimated using the 
Z-Peg code by integrating the total star formation his- 
tory (SFH) of the best-fit model, up to the best-fit age 
and subtracting off mass loss from late stages of stellar 
evolution. 

3. CLASSIFICATION OF WET AND DRY MERGERS 

We used two techniques to try to distinguish between 
wet and dry mergers. The common starting point for 
both methods is subdivision into red and blue stellar 
populations on the basis of rest-frame color relative to 
a fiducial reference color (Van Dokkum 2005; Willmer et 
al. 2006). Our approach to defining this reference color 
uses the rest frame g' — i' versus g' color-magnitude di- 
agram for ~ 2200 visually classified field galaxies in the 
CFHTLS Dl and D2 field (see Figure [1]). Red dots rep- 
resent visually classified elliptical galaxies and blue dots 
indicate spiral galaxies. The green dots indicate the ob- 
jects with MIPS 2Anm detection (down to a fiux limit of 
340 ^Jy). iCowie fc Bargei] ()2008f ) report that at z < 1.5 
most red galaxies with a 24/xto flux > 80/iJy fall into the 
blue cloud after the appropriate dust extinction is ap- 
plied. To account for the color change in dusty sources, 
we artificially assign the green dots with g' — i' color 
greater than the color bimodality to the blue cloud. In 
addition, color bimodality is known to be magnitude or 
stellar mass-dependent and us ually derived from the fit- 
ting of red sequenc e objects (|van Dokkum et all [2000l: 
I Willmer et al.l [20061 ) . Therefore, the following stellar- 
mass-dependent fiducial color cut was adopted based on 
the red sequence fitting of 115 visually classified ellipti- 
cal field galaxies with spectroscopic redsh ifts from Cos- 
mological Evolution S urvey (COSMOS) ([Scoville et al.l 
120071: iLillv et al.ll2007[) . The fitting line is expressed by 
the following equation: 



= -0.0076 + 0.13 X Mstellar ^ C (!) 



The constant C serves as a parameter to control the verti- 
cal position of the fitting line on the diagram. To account 
for the potential classification errors caused by different 
slopes and C values, we have explored the implications 
of changing the free parameters in Equation (1) and find 
that all trends reported in this paper remain robust to 
the specific numerical values chosen. To show this, in 
several key figures in this paper we will bracket the re- 
sults obtained using Equation (1) with curves showing 
the envelope obtained by using three values of C (0, 0.1 
and 0.2) as well as using a perfectly horizontal line set at 
g' — i' = 1.29 mag to divide field galaxies into blue and 
red cloudt|3 

3.1. Method I: integrated colors 

The first method for segregating wet mergers from dry 
mergers is, essentially, the simplest conceivable. We look 

^ The reason to use g' — i' = 1.29 is that we try to balance the 
contamination of blue (red) galaxies in red (blue) cloud. For the 
g' — i' = 1.29 color cut the contamination on both sides is 14%. 
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at the total integrated color of the merger which we 
obtained by summing over all pixels in the merger (or 
merging pair), and simply note if the integrated color 
of the complete system is redder or bluer than the fidu- 
cial color threshold. Mergers whose integrated colors are 
redder than the threshold are deemed 'dry', and systems 
bluer than the threshold are deemed 'wet'. As mentioned 
above, the only caveat is that we re-assign dry mergers 
with 24/^m detections to the wet merger category on the 
assumption that these are blue galaxies being reddened 
by dust. 

This simple method is straightforward but it is not at 
all obvious that loss of the spatial information is an ac- 
ceptable trade-off for such simplicity. For example, what 
if one object in a merging system is redder than the fidu- 
cial color, while the other object is bluer? We therefore 
decided to explore a somewhat more refined approach 
to 'wet' vs. 'dry' merger classification that retains some 
component of the spatial information in the images. 

3.2. Method II: spatially resolved colors 

Our second method is based on analysis of the colors 
of individual pixels. Pixels with rest-frame g' - i' color 
greater than the fiducial threshold are labeled as 'red', 
and the ratio of the total fiux in red pixels to the fiux in 
all pixels is calculated. We refer to this quantity as the 
'Red-to- Total ratio', (R/T), given by: 

{R/T) - ^ (2) 

^ total 

where Fred indicates fiux contained in red pixels and 
Ftotai refers to the flux from the entire merger. 

The idea behind (R/T) is to crudely decompose a 
merger into individual stellar populations. A galaxy with 
[R/T) = 1 is comprised entirely of stellar populations 
with colors on the red sequence, while a merger with 
{R/T) = is comprised entirely of stellar populations in 
the blue cloud. Of course most mergers are expected to 
lie somewhere in between these extremes. For the sake of 
simplicity we choose to define systems with (R/T) > 0.5 
as being 'dry', but throughout this paper we will explore 
the results obtained for range of {R/T) values to ensure 
that our conclusions are not tied to any specific value of 
{R/T). Finally, we note that in common with Method I, 
we re-assign 'dry' mergers with MIPS 24 /im detections 
into the blue cloud. 

4. RESULTS 

4.1. Merger Fractions 

Figure [3] shows the merger fraction for a stellar mass- 
hmited sample of objects (M* > 10^-^ Mq), /m, as a 
function of redshift. We define this merger fraction in the 
obvious way as /a/ = Nm/N, where Nm is the number of 
mergers and N is the number of galaxies in the complete 
sample (once again, we refer the reader to Bridge et al. 
2010 for details of the full sample). The figure shows 
curves for wet mergers (blue lines), dry mergers (red 
lines), and all mergers (black lines). The left-hand panel 
shows the merger fraction computed using integrated col- 
ors while the right-hand panel shows the corresponding 
curves computed using resolved colors, at three differ- 
ent values of the red-to-total fraction used to segregate 



wet and dry mergers. Errors were estimated assuming 
Poisson statistics. 

A number of interesting results emerge from these fig- 
ures. Firstly, a comparison of the left and right panels 
shows that the general trends seen are quite indepen- 
dent of the specific methodology used. In all cases the 
dry merger fraction is almost equal to the wet merger 
fraction at the lowest-redshift bin (around z ~ 0.2) but 
the fraction of dry mergers decreases rather quickly be- 
tween z = 0.2 — 0.3 and remains somewhat flat at higher 
redshifts. In contrast to this, the fraction of wet mergers 
increases rather rapidly with redshift, and by z 0.7 wet 
mergers outnumber dry ones by a factor of 6 to 1. The 
dry merger fraction atz = 0.5is^l%, which is in agree- 
ment with t he dry merger fraction of 1% at (z) ~ 0.55 
obtained bv lDe Propris et al.l ()2010D using close pairs. 

Another interesting feature noted in Figure[3]is the way 
that the total merger fraction drops from around 6.5% 
at z ~ 0.2 to around 4.5% between z = 0.3 and z = 0.5, 
before rising to around 7% at z = 0.7. This effect seems 
due to the rapid growth with time in dry mergers at 
lower redshifts and the rapid decline with time in wet 
mergers in earlier Universe. This dip is seen more pre- 
dominantly in D2 than i n Dl field sugges ting that cos- 
mic variance is involved (|Bridge et al.li201Q) . The non- 
monotonic changes in the merging fraction seen suggest 
that over the redshift range being probed a power law 
provides a poor description for the changing total merg- 
ing fraction. At the same time, the monotonic changes 
in the wet and dry populations when considered sepa- 
rately suggests that a power law might provide a reason- 
able description for these sub-populations, so we have 
analyzed them separately. We adopt the analytic form 
/m = /o X (1 + z)"^, where /o is also set as a free parame- 
ter, and m is the standard power-law index. Fitting this 
relationship to the R/T = 0.5 curves on the right-hand 
panel yields m = 2.0 ± 0.3 and m = —3.1 ± 0.5 for wet 
and dry mergers, respectively. Since the merger sample 
may not be complete in the last bin at z = 0.77, and 
the power-law indices m = 1.2 ± 0.4 and m = —4.1 ± 0.6 
are obtained for wet and dry mergers with the last data 
point excluded. We see that the dry merging population 
is best described by negative evolution, and the basic 
statement that the wet and dry populations are evolving 
with the same m is excluded at more than three sigma 
significanceQ. 

The careful consideration given to selection effects de- 
scribed in Bridge et al. (2010) leaves little doubt that 
the effects noted in Figure [3] are really seen in the obser- 
vational data, but what is their ultimate physical signifi- 
cance? The visibility of red and blue stellar populations 
is expected to be a strong function of k-corrections, which 
differer markedly for red and blue galaxies at z ~ 0.7 
(the redshift limit for the present paper). To better un- 
derstand this, we will now explore the space densities of 
mergers using the V^nax formalism to account for color 
and luminosity biases. 

4.2. The Space Density of Merging Galaxies 

Note that the total and wet merger fraction keep increasing 
rapidly beyond the redshift limit of this paper, so it seems that the 
merger rate evolves more dramatically at higher redshifts than at 
z < 0.7. See Bridge et al. (2010) for more discussion on this point. 
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Fig. 3. — Merger fractions computed using two different methods. The left panel shows the merger fraction derived from different fiducial 
color cuts with integrated colors, while the right panel shows the merger fraction derived from different internal color ratios (see text for 
details). It is obvious that different color cuts do not affect the scientific results. In both panels, red curves indicate dry mergers, blue 
curves indicate wet mergers, and black curves indicate the total merging population. It is clear that the increase with redshift in the wet 
merger fraction, and the decrease with redshift in the dry merger fraction, are robust and do not depend on the selection method. Error 
bars are estimated by assuming Poisson errors. 




Fig. 4. — Space densities of wet and dry mergers. As in the previous figure, blue and red curves correspond to wet and dry mergers, 
respectively, while the black curve shows the total for all mergers. The growth in the total space density of mergers is modest, increasing 
by about a factor of two over the redshift range probed. The wet and dry merging galaxies show opposite trends with redshift. The space 
density of wet mergers is increasing with redshift, while that of dry mergers is (perhaps) modestly decreasing. 



Figure S] shows the evolving space density of all merg- 
ers, as well as of the sub-samples of wet and dry merg- 
ers. These space densities were computed by weighting 
each merger by 1 / Vmax and summing ov er redshift bins in 
the st andard manner (see, for example, iTotani fc Yoshiil 
(|1998D ). As described in §2.2, V;„ax values were com- 
puted using the Z-PEG code, although it is once again 
emphasized that in most cases the ultimate limit on the 
detectable volume is set by the merger classification limit 
rather than the magnitude detection limit. As with the 
previous figure, the blue and red curve correspond to wet 
and dry mergers, while the total for the whole popula- 
tion is shown in black. Error bars were computed by 
bootstrap resampling. 

The main trends noted in our description of Figure |3] 
remain visible in Figure |4j most notably the mild nega- 



tive evolution in the dry mergers and positive evolution 
in the wet mergers, with the populations converging at 
around z ~ 0.2. However, a rather striking new feature 
emerges: the total space density for all mergers appar- 
ently remains almost flat at all redshifts, with a mild in- 
crement at z > 0.5. In other words, the decline with 
redshift in dry mergers appears almost perfectly offset 
by a rise with redshift in the fraction of wet mergers. 

4.3. The Role of Stellar Mass 

As was noted in the Introduction, stellar mass is 
the central parameter driving much of galaxy evolution. 
What role does stellar mass play in conditioning the ef- 
fects noted in the previous two figures? To explore this, 
we applied a stellar mass cut of 10^'^Mq and a redshift 
cut of z > 0.1 to our merger sample, leaving 1296 ob- 
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Fig. 5. — [Left:] Mass density function in mergers in two redshift bins. The high-redshift bin (0.4 < z < 0.7) is shown with a dashed 
curve, while the low redshift bin (0.1 < 2; < 0.4) is shown with a solid curve. We see evidence for modest evolution in the mass density 
function, with most change occurring at intermediate masses, and no change at the high and low mass ends. [Right:] Subdivision of the 
data in the left hand panel into wet and dry mergers. Line styles segregate the data into two redshift bins, as for the left-hand panel. As 
in previous figures, red curves correspond to dry mergers, and blue curves correspond to wet mergers. For mergers in the most massive 
bin, an increase in the dry merger space density is offset by a decrease in the wet merger space density, so the total space density is nearly 
conserved. At intermediate masses, the mass density function of dry mergers is nearly unchanged in both redshift bins, with perhaps some 
evidence for a slight decrease in the space density of intermediate-mass dry mergers at high redshifts. On the other hand, the mass density 
function of wet mergers is increasing with redshift. 



jects in our mass-limited sample. Figure [5] shows the 
stellar mass density function for this total sample in two 
redshift bins (left-hand panel), and the result obtained 
when the total sample is divided into wet and dry sub- 
populations (right-hand panel; note we have adopted the 
same color scheme as for the previous figures). Error 
bars were derived in the same manner as for the merger 
number density. 

The left-hand panel of Figure [5] shows that the in- 
crease in the space density of high-redshift mergers oc- 
curs mostly in intermediate-mass galaxies. There is little 
evidence for an increase in the space density of merging 
galaxies in either our lowest mass bin or in our highest 
mass bin, which is consistent with the result shown in 
Figure |4l The overall situation is clarified further in the 
right hand panel, which shows that: (i) the space density 
of wet mergers is always higher than the space density of 
dry mergers, regardless of redshift or mass; (ii) the space 
density of wet mergers evolves more quickly than the 
space density of dry mergers. In fact, in the dry merging 
population there appears to be little evidence for redshift 
evolution except in the lowest mass bin, suggesting that 
any evolutionary trends for dry mergers occur mainly at 
low masses. 

5. DISCUSSION 

We would like to compare results with those presented 
in recent papers by Lin et al. (2008), hereafter Lin08, and 
Jogee et al. (2009), hereafter Jog09. These authors used 
different selection methods. LinOS used close pair selec- 
tion while JogOQ used morphologi cal criteria quantified 
using the automated CAS method (jConselice et al.|[2000l : 
iConse licc 2003). Both sets of authors provide estimates 
of the volumetric merger rate with respect to redshift, 
which provides a nice focus for inter-comparisons. The 
Lin08 merger sam ple is composed of data from DEEP2 
(iDavis et al.ll2003D. the T eam Keck Treasury Redshift 
Survev (^^yirth et al.ll2004D. th e Southern Sky Redshift 
Survey 2 Cda Costa et aL'^QDS*), the IVIill ennium Galaxy 
Catalog ((Liske et al., ,2003; ,Driver et al.l 12005. ) and the 



Canadian Network for Obs ervational Cosm ology Field 
Galaxy Redshift Survey 2 (|Yee et al.l[2000l ). for a to- 
tal of 506 close pairs (compared to 1300 mergers pre- 
sented in thi s work). The J og09 sample is taken from 
the GEMS ()Rix et al.l I2004D with spectrophotometric 
redshifts and spectral energy distributions taken from 
COlVIBO-17 (I Wolf et al.ll2Q0l . for a total of ~800 merg- 
ers. 

In order to compare our results with Lin08 and Jog09, 
the merger densities shown in Figure|4]must be converted 
to a merger rate density, and this requires a merging 
timescale to be assumed. The appropriate timescale to 
use is unclear, since it may be a function of a number 
of factors, such a s galaxy mass ratio and gas fraction 
(jLotz et al.ll2009bl ial). Additionally, different merger se- 
lection methods would select meregrs in different merging 
stages and therefore require different correction factors 
as well as merging timescales. To make sure that we 
are comparing the same thing in Figure |6l we still keep 
the original merging timescales used in the corresponding 
papers. Recent work based on N-body/hydrodynamical 
simulations for equal mass gas-rich mergers suggests that 
merging timescales typically range from 0.2 Gyr to 0.9 
Gyr, on the basis of a comparison between the tidal fea- 
tures seen and those produced by simulation^ Note 
that the simulation results are valid for mergers selected 
by different methods. At present we will adopt a merging 
timescale of 0.8±0.2 Gyr estimated from visual compar- 
ison of both dry and wet mergers presented in Bridge et 
al. (2010) who also visually examined the snapshots of 
simulated mergers and the duration in which the galax- 
ies would be classified as interacting' under the criteria 
used to classify the galaxies presented in the paper. Fig- 
ure|6]shows the comoving merger rate density for galaxies 
with redshifts z < 0.7 in these samplefl In this figure 

^ It is worth noting that the most obvious features emerg e at the 
first encounter and the final merging stages (Lotz et al. 20 08al ) . 

* LinOS designate around 30% of their merger sample as 'mixed 
pairs' (one red galaxy and one blue galaxy). For purposes of com- 
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Fig. 6. — A comparison of the evolving merger rate density in the present work with the results obtained by Lin et al. (2008; LinOS) 
and Jogee et al. (2009; Jog09) using pair counts and automatic morphological selection. Our merger rate density is derived by choosing 
a merging time scale of 0.8 ±0.2 Gyr, as discussed in Bridge et al. (2010). The merging timescale used by Lin08 and Jog09 is 0.5 Gyr. 
Qualitative trends determined using both methods are similar. The absolute merger rate density from our work is lower by a factor of 
compared to LinOB but is comparable to that in Jog09 (see text for details). 



dashed lines show the resuhs from LinOS, the dotted hne 
corresponds to Jog09, and sohd hnes show results ob- 
tained from our data. Blue and red colors represent wet 
and dry mergers, respectively. Note that Jog09 does not 
break the sample down into wet and dry mergers, so only 
the total merger rate density is shown. 

Figure [6] shows several interesting things. Firstly, we 
see that the merger rate density from our work is gener- 
ally a factor of ~ 3 less than that in LinOS, although it is 
in excellent agreement with that of Jog09. Note however, 
that our merger fraction (shown in Figure [7]) is actually 
in rather good agreement with that of LinOS, so that 
the disagreement between our merger rate density and 
that of LinOS is due to differences in translating from an 
empirically close pair fraction to a merger rate density. 
This step is more complicated in a pair count analysis 
than it is in a morphological analysis, because the for- 
mer involves a number of coefficients (fraction of galaxy 
pairs to mergers and correction factor for the compan- 
ion selection effect due to the limited luminosity range) 
with values that are not well established. The relevant 
issues are well-described in LinOS, who outline the steps 
needed to obtain Equation 5 in that paper. Another fac- 
tor that may explain the difference with regard to LinOS 
include the possibility that we are missing some mergers 
because their tidal tails have dropped below our surface 
brightness detection threshold, although this possibility 
has been investigated in Bridge et al. (2010) and we be- 
lieve it is not a significant factor (i.e. certainly less than 
a factor of two) out to z = 0.7. 

Another interesting result that emerges from the fig- 
ure is the fact that in spite of the systematic offset in 

parison with our own data, we have simply included half of the 
LinOB galaxies in the wet mergers category and half of them in the 
dry mergers category. None of the conclusions in this paper are 
sensitive to this assumption. 



the absolute rates compared to LinOS, many of the qual- 
itative trends are similar when our results are compared 
to theirs. In both cases, the density of wet mergers is 
essentially identical to the density of dry mergers in the 
nearby Universe. Also in both cases the evolution in the 
total merger rate densities at these redshifts is at best 
modest, and arguably flat. Finally, in both cases we see 
that wet mergers and dry mergers have different trends 
as a function of redshift, and that dry mergers contribute 
less at higher redshifts than they do at lower redshifts. 
Since Jog09 do not break their sample into wet and dry 
mergers, we are unable to compare such trends with their 
sample, although as we have already noted our overall 
merger rate densities are in good agreement with theirs. 

It is highly interesting to explore how the trends shown 
in Figure [6] depend on the stellar masses of the merg- 
ers. This can be most clearly seen by placing the stellar 
mass density functions for wet and dry mergers (already 
shown in Figure [5l subdivided into two redshift bins) 
into a broader context, namely that which includes the 
population these systems are thought to develop into in 
hierarchical models. To this end. Figure [8] shows the 
mass density function (converted to our chosen IMF) for 
early-type galaxies at both z ~ 0.4 and z ~ 0.7 taken 
from Bundy et al. (200^, as well as our derived mass 
density functions for mergers. If early- type galaxies are 
being built up via mergers then the mass density func- 
tion of early type galaxies is a cumulative quantity that 
grows with redshift. A simple model that is worth con- 
sidering in one in which any redshift-dependent trends 
have no mass dependence. On a logarithmic plot such a 
uniform multiplicative mass growth simply displaces the 
mass density function upward, so the normalization of 
the mass density function changes with redshift, but the 
shape of the curve stays the same. Although this argu- 
ment is based on an assumption that the merger products 
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Fig. 7. — Comparison of merger fractions witii two different merger selection metliods. Tlie blue line and triangles indicate the merger 
fraction derived from close pair merger selection method; the red line and circles indicate the merger fraction derived from mergers selected 
by automated CAS asymmetry and dumpiness parameters. In Figure |6] our merger rate density is in good agreement with data obtained 
by Jog09. However, it is clear that our merger fraction is in good agreement with data obtained by LinOS. This indicates that the main 
uncertainty comes from different merging timescales and correction factors that translate the empirical merger fraction to merger rate 
density (see text for more details). 

(~ IO^^'^Mq), which outnumber mergers by a factor of 
about 50 in Figure [H In this case, a merging timescale 
of order 0.5 Gyr is consistent with mergers producing 
the space density of ehiptical galaxies seen in Figure [5] 
(again, assuming a constant merger rate, as is suggested 
from our data). 

We conclude from this analysis that not all massive 
major mergers (Mstciiar > lO^^Mo) will end up being the 
most massive early-type galaxies, at least not at z < 0.7. 
On the other hand, the merging timescales suggested by 
numerical simulations are consistent with the data if we 
suppose that less massive early-type galaxies form via 
mergers. Since low-intermediate mass ellipticals are 10- 
100 times more common than their most massive coun- 
terparts, the hierarchical explanation for the origin of 
early-type galaxies may be correct for the vast major- 
ity of early-types, even if incorrect for the most massive 
ones. 

If the massive merging systems captured in Figure [8] do 
not end up as early- type galaxies, what is their ultimate 
fate? That is unclear. Howe v er. [Robertson k, Bullock! 
(|2008D : [G"overnato et al.l (|2009[ ): lKannappan et al.l(|2009( ) 
claim that the end product of major mergers of gas-rich 
spiral galaxies is often a larger spiral galaxy, so that 
would seem to be a strong possibility. 

6. CONCLUSION 

We have divided the 1298 merging galaxies in Bridge 
et al. (2010) into 'wet' and 'dry' merging categories us- 
ing two techniques, one based on integrated colors and 
the other based on spatially resolved colors. We find 
that our general results are independent of the specific 
method used. Using the VJnax methodology, we are able 
to compute the space density of merging galaxies out to 
redshift z ^ 0.7. The local space density of wet mergers 
is essentially identical to the local space density of dry 



do not participate in subsequent mergers, since the typ- 
ical merger rate is ~ 10~^ Gyr^^Mpc"'^ it is not very 
likely that a typical early-type galaxy will experience 
multiple massive mergers within the span of a few bil- 
lion years. It is clear from Figure|8]that the shapes of the 
mass density function curves are completely different and 
that of early type galaxies has very little evolution at the 
high mass end (with observational errors). This seems 
rather surprising in light of the expectation that major 
mergers are thought to be the progenitor population for 
early-type galaxies. In that case, one expects that the 
two curves would have similar shapes, with the early- 
type galaxy mass density function displaced upwards as 
elliptical galaxies form, unless the merging timescale is 
a very strong function of mass. We conclude from this 
analysis that if early- type galaxies are formed from merg- 
ers, at some epoch the mass density function of mergers 
must have been very different from that seen today. 

We can explore the implications of the very different 
forms for the merger and early-type galaxy mas s density 
functi ons by adopting the approach first used by iToomrd 
(|1977D , who compared the number of early- type galaxies 
to the numbe r of merging galaxies in the Second Refer- 
ence Catalog (|de Vaucouleurs et al.lll976l ). At the high- 
mass end, we find that the space density of massive early- 
type galaxies is only a factor of five larger than the space 
density of wet mergers. This is a highly puzzling result to 
understand in a framework in which massive mergers pro- 
duce elliptical galaxies, since our results also show that 
the merger rate has not changed rapidly since z=0.7. Un- 
less the merging timescale is ~ 5 Gyr for massive merg- 
ers (about an order of magnitude larger than current 
estimates) the space density of massive galaxy mergers 
would greatly over-produce massive elliptical galaxies if 
these mergers all turn into early-type systems. A better 
match is seen for elliptical galaxies of intermediate mass 
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Fig. 8. — A comparison of the stellar mass density function for mergers in the present sample with the corresponding stellar mass density 
function for early- type galaxies presented in Bundy et al. (2005). The data point of early type galaxies are more or less from the same 
redshift bins. Two vertical lines indicate the comp leteness l i mit o f early type galaxies in different redshifts. The x-axis is in units of 
Mstcllar h^/Msun to facilitate comparison with Bundy et al.l II2005I V Note the completely different shapes for these two samples, even 
though elliptical galaxies are associated with the end product of mergers, and there is no evolution in the space density of massive elliptical 
galaxies. This shows quite clearly that not all of the major mergers will end up with massive early-type galaxies, unless the merger timescale 
is much longer than usually assumed. See text for details. 



mergers. The space density of merging galaxies does not 
change rapidly with redshift, increasing by less than a 
factor of two out to z ~ 0.7. On the other hand, wet 
and dry merging populations show different evolutionary 
trends. At higher redshifts dry mergers are a smaller 
fraction of the total merging galaxy population, while 
the wet mergers make a proportionately greater contri- 
bution. 

We have compared the stellar mass density function of 
merging galaxies to the corresponding stellar mass den- 
sity function of early- type galaxies, and find that these 
functions are very different in shape, even though the for- 
mer are thought to be the progenitor population of the 
latter. We also show that the present space density of 



massive galaxy mergers is already similar (within a fac- 
tor of three) to that of existing massive elliptical galaxies. 
This suggests that not all of the massive major mergers 
will end up with massive early-type galaxies, unless the 
merging timescale is much longer than expected. On the 
other hand, for systems with masses less than IO^^'^Mq, 
we find that the space density of low-intermediate mass 
elliptical galaxies is consistent with hierarchical forma- 
tion via mergers. 
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